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Abstract. The objective of this paper is to present a modern and concise new deriva¬ 
tion for the explicit expression of the interior and exterior Newtonian potential generated 
by homogeneous ellipsoidal domains in (with A ^ 3). The very short argument is 
essentially based on the application of Reynolds transport theorem in connection with 
Green-Stokes integral representation formula for smooth functions on bounded domains 
of R^, which permits to reduce the A-dimensional problem to a 1-dimensional one. Due 
to its physical relevance, a separate section is devoted to the the derivation of the demag¬ 
netizing factors of the general ellipsoid which are one of the most fundamental quantities 
in ferromagnetism. 


1 Historical Introduction, Motivations 

The computation of the gravitational potential induced by an homogeneous ellipsoid was one 
of the most important problems in mathematics for more than two centuries after Newton 
enunciated his universal law of gravitation n El El IMl |3Q1 E]. 

Once fixed an ellipsoidal domain D of IR^, the problem consists in hnding an explicit 
expression for the Newtonian potential induced by a constant mass/charge density on 

M- 

in the internal points of {interior problem) and in the exterior points of {exterior prob¬ 
lem). 

In the case of a homogeneous spherically symmetric region, Newton (in 1687) proved 
what nowadays is known as Newton’s shell theorem 131121]: if n is an homogeneou^ 
spherical region centered at the origin, then for all t > 1, is constant in fl, 

i.e. tfl\fl (the so called hollow ball) induces no gravitational force inside Q. Moreover, a 
spherically symmetric body affects external objects gravitationally as though all of its mass 
were concentrated at a point at its center. 

For what concerns an ellipsoidal domain of the ellipsoidal interior problem was for 
the hrst time solved by Gauss (in 1813) by the means of what it is in present days known 
as the Gauss divergence theorem |16) . Later, in 1839, Dirichlet proposed a solution of 
the interior problem based on the theory of FOURIER integrals m- 

The results of Gauss and Dirichlet can be summarized by saying that if is an 
ellipsoidal domain centered at the origin, then the gravitational potential induced by an 
homogeneous ellipsoid in its internal points is a second order polynomial. In other terms: 

Mn[^Q\{x) = c - Px ■ X Vx E n, (2) 

^By an homogeneous region of we mean an open and connected subset of endowed by a constant 
density of masses/charges. 




for some constant c € IR and some matrix P G whose values can be expressed in terms 
of elliptic integrals jlHl |35] . 

Remark 1.1. For the sake of completeness we recall that the converse statement (the inverse 
homogeneous ellipsoid problem) is also true |11[ ITSl [T8] . namely: if 17 is a bounded domain of 

such that IR^\n is connected and ([2|) holds, then Q is an ellipsoid. Historically speaking 
the inverse homogeneous ellipsoid problem was for the first time solved by DiVE |12) in 1931 
for N = 3 and in 1932 by HOLDER m for N = 2. A modern proof of this result can be 
found in DiBenedetto and Friedman m who, in 1985, extended it to all N ^ 2. In 
1994, Karp |18) . by the means of certain topological methods, obtained an alternative proof 
of the inverse homogeneous ellipsoid problem. 

Despite the well-knownness of (|2]) in the mathematical and physical community, and its 
importance in theoretical and applied studies |Il[5l|7l[32l|33llinil25], rigorous proofs of that 
result are not readily available in the literature: to the best knowledge of the author, relative 
modern treatments of the interior problem can be found in m and |35) . and more recently 
in |22) and |26) where also the the exterior problem is investigated. However, in all the cited 
references, the solution of the problem is always based on the use of ellipsoidal coordinates 
which tends to focus the attention on the computational details of the question rather than 
on its geometric counterpart. Eventually, modern proofs of Newton’s theorem and relation 
([2]), as well as far reaching beautiful generalizations, can be found in [19] where neverthless 
the problem of finding an anlytic expression for the coeficients c and P is not touched. 

Aim of this paper is to give a modern and concise derivation for the expression of the 
interior and exterior Newtonian potential (induced by homogeneous ellipsoids). The very 
short argument is essentially based on the application of REYNOLDS transport theorem in 
connection with Green-Stokes integral representation formula for smooth functions on 
bounded domains of R'^. This approach permits to reduce the A^-dimensional problem to a 
1-dimensional one, providing (in particular) at once a proof of ([2|) together with an explicit 
expression of the coefficients P and c in terms of 1-dimensional integrals. More precisely, 
the paper is organized as follows: 

Section [2] is devoted to the main result of the paper. We give a concise proof of the homo¬ 
geneous ellipsoid problem. For completensee, we then derive Newton’s shell theorem 
as an immediate corollary. 

In Sectior(3]we focus the attention to the three dimensional case. An expression in terms 
of the elliptic integrals is given for the coefficients of P. Due to its physical relevance, par¬ 
ticular attention is paid to the eigenvalues of P. Indeed, when N = 3, the matrix P and 
its eigenvalues, known in the theory of ferromagnetism respectively as the demagnetization 
tensor and the demagnetizing factors, are one of the most important and well-studied quan¬ 
tities of ferromagnetism n El [510 dm [25]. In fact, the following magnetostatic counterpart 
of the homogeneous ellipsoid problem holds: given a uniformly magnetized ellipsoid, the in¬ 
duced magnetic field is also uniform inside the ellipsoid. This result was for the first time 
showed by POISSON |28j . while an explicit expression for the demagnetizing factors was for 
the first time obtained by MAXWELL |23j . Their importance is in that they encapsulate the 
self-interaction of magnetized bodies: their knowledge being equivalent to the one of the 
corresponding demagnetizing (stray) fields |5]- 

2 The Interior and Exterior Potential of an Homogeneous Ellipsoid 

In what follows we denote by H the ellipsoidal domain of R'^ {N ^ 3) having ai, 02 ,..., oat 
as semi-axes lengths. We then denote by {flt)te[o,+oo) the family of ellipsoidal domains of 


, given by the inverse image ^{B^) of the unit ball of IR^ under the one parameter 
family of diffeomorphisms (j)t : x £ i—?■ y/~^x G IR^, where 


:= diag 


V^ofTl sjal + t + t 


(3) 


Note that each diffeomorphism ^ maps the unit ball of onto the ellipsoidal domain 
of R^ dehned by the position fit = {x G R^ ; \<j)t{x)\'^ ^ 1}. In particular dQ.t = {x G 
R^ : \M x)p = 1} and fl = fig. Finally, we denote by A/nJlnJ the Newtonian potential 
generated by the uniform space density of masses or charges on 

MnA^nt]{x) = CN j I-^7^^dr(?/), (4) 

Jnt I® “ y\ 

with cjv:=[(-N — 2)ccAr]“^ and the surface measure of the unit sphere in R-^ (cfr. |13[I14] 1. 
The main result of the paper is stated in the following: 

Theorem 2.1. Let Q = {x G R^ : |(/>o(x)p ^ 1} he the ellipsoidal domain ofR^ having 
(ai, 02 ,, Oat) G as semi-axes lengths. For every x G R^ 


r+oc 

A/o[lo](x) = - / -ft{l - Atx ■ x)dt, 

^ Jt{x) 


N 


:=n 

2=1 


Qi 


where we have denoted by r the non-negative real valued function 


(5) 


f 0 if X £FL 

\ r(x) if X G n 


( 6 ) 


Remark 2.1. For general ellipsoidal domains the integral in ([5|) can be evaluated by the 
means of the theory of elliptic integrals (cfr. |T1129]). The function r dehned by ([6| can be 
computed by solving the equation = {x G R'^ : |(/)T-(a,)(x)p = 1}. In particular, 

when n is a spherical region of radius a (centered around the origin), the function r reduces 
to the function equal to |xp — o^ if x G R^\n, zero otherwise, and the integral in (j5]) can be 
readily computed. 


Proof. For every t G R"*", the function — 1 is the unique solution of the homogeneous 
Dirichlet problem for the POISSON equation u = 2tr(^t) in Lit, u = 0 on dLlf Thus, 
according to the Green-STORES representation formula (see mi), for every t G R"*" we have 

(|(/)tp - l)lf2t = -2tv{At)J\fn,Ant] + SanAdn\4>t\A “ R^\dLlt, (7) 


where we have denoted by 

<Sdnt[dn.\f>tWx) := CN [ ^ da{y), (8) 

Jaut \x-yr 

the simple layer potential generated by the space density of masses or charges dn\4>t\‘^ con¬ 
centrated on dLlt (cfr. |13[ 114) 1. Next, we observe that the A^-dimensional Newtonian kernel 
Ixp”'^ is in kFjoc (R'^)) therefore due to Reynolds transport theorem, 

=CN [ Z^-2 My) = \sanAdM% 

Jant ~ y\ ^ 


( 9 ) 















where we have denoted by Vt := dt[(j)i o (fi^ the eulerian velocity field associated to the 
motion for which one has 

vt{y) ■= dt[(j)i^]o (l)t{y) = ^Aty = ^V\4>t{y)f V(t,y) E [R+X (9^4. (10) 

Hence, substituting ([9]) into ([7} we get 

\{\Mx)\^ - l)lnt(a;) = -^tT{At)J^n,[lnt]{x) + dtAfnt[lnt]ix) Vx E R^XdQf (H) 
Moreover, by the assignment 

7 t:=exp(--/ tr(H^)ds)=JJ ^ 7o = 1, (12) 

V 2 Jo J 

the equality (II ip reads as 

dt{-itNnt[^nt]) = VxE[R^\(9Ht. (13) 

Thus, once introduced the non-negative real function defined by ([6| we have, for every t E R"*" 
and for every x E R^ 

dt{ntNQ.,[lnt]){x) = - l)l[r(x),+oo)(t)- (14) 

Integrating both members of (1141) on [0, -|-oo); taking into account that due to the well-known 
decay at infinity of the Newtonian potential |34) one has limt_).+c« yjA/bJln^] = 0; we finish 
with (j5|) . □ 

Corollary 2.1. (Newton’s shell theorem) Let Q C be an homogeneous spherical 
region (centered around the origin) of radius a and of total mass M. For every x E R^\H the 
induced gravitational potential is the same as though all of its mass were concentrated at a 
point at its center. Moreover, for all t > 1, A/'tr2\n[ltr2\n] constant in Ll, i.e. the hollow ball 
induces no gravitational force inside Ll. 

Proof. We denote by p1q{x), {p := M/|H|) the uniform density of mass in H. The gravita¬ 
tional potential induced by p in R^ is then given by Up = dvrGA/n)/?] = dTrGA/nflo], where 
we have denoted by G the gravitational constant. In this geometrical setting the function r 
defined by ([6| reduces to the function r(x) := (|xp — a^)lRiv\n(x) and the integral in ([5)) 
immediately gives 

Up{x) = -ttGp ( -if X E H , Up{x) = -vra^-—^ = -—— if x E R^\H. (15) 

3 \ 2 / 3 |x| |x| 

Thus, for every x E R^\H the induced gravitational potential is equal to the one induced 
by a Dirac mass concentrated in the center of D. The gravitational field is given by g := 
—'Vup and the fact that the hollow ball tLl\Ll, t > 1, induces no gravitational force inside 
D can be immediately seen by splitting the uniform density of mass in tLl\Ll in the form 
p = {M/\tLl\Ll\)ltn — (M/|tD\D|)lQ. Indeed by linearity we get that Up is constant in D 
and therefore g = 0. □ 





3 The demagnetizing factors of the general ellipsoid 


We know focus on the three-dimensional framework (iV = 3) and, in particular, on the 
so-called demagnetizing factors of the general ellipsoid |25) . To this end we recall that the 
demagnetizing (stray) field h\m] associated to a magnetization m G C°°{ri) can be expressed 
as the gradient field of a suitable magnetostatic potential (p^n (see [ISIET]). Precisely tprn '■= 
— divA/n[iLi] and h\rn\ := —Vipm in In particular, if m is constant in fl, then = 
—m • VA/h)!^]. Thus from ([5]): 

(Pm{x) = Px ■ m , h[m\ =—Pm Vx G fl, ( 16 ) 


where we have denote by P := VA/b[lQ] the diagonal matrix, known in literature as the 
demagnetizing tensor, whose diagonal i-entry (the i-th demagnetizing factor), by virtue of 
m, is given by 


Pi:= 


1 

2 




Vi G IM3. 


(17) 


Proposition 3.1. We have Pi ^ 0 for every i G IKI3 and if ai ^ 02 ^ 03 then Pi ^ P 2 ^ P3. 
The trace of P satisfy the relation tr(P) = 1. 


Proof. The first statement is obvious. The relation tr(P) = 1 can of course be verified 
by a direct evaluation of the integrals in (I17p . but it is also possible to observe that since 
the Newtonian potential satisfies the PoiSSON equation AA/b[lfi] = — lo, one has 

tr(P) = div(Px) = -AA/b[lf 2 ] = 1^. □ 

Assuming ai ^ 02 ^ 03 , from m and the theory of elliptic integrals we get 


Pi = 1 — P 2 — P 3 
as 


P2 = — 

( 

P 3 = +- 


I -«i 

a 2 


as 


0102 


Oo o 


3 L 


0-2 


{aj - {; 

O 1 O 3 


02 

arccos | — 

Ol 


of -a| 


ot — a 


2 / J 


(ai-ai)V2-(^ 


TP I I 

P I arccos | — 

Ol 


of - oi 


ot — o 


s / J 


where, for every y G K and every 0 < p < 1 we have denoted by 


P(y\p} ■= [ (1 — psin^ 0)^/^d0, 

Jo 


(18) 

(19) 

( 20 ) 


( 21 ) 


the incomplete elliptic integral of the second kind expressed in parameter form II [29]. In 
particular, in the case of a prolate spheroid (oi ^ 02 = Os) we get 


Pi = - 


(of - [ 

l-Fi 


(of — of )^^^ -|- Ol arccoth 


Ol 


(of - oi)V 2 


P 2 = P3 = 

z 

while in the case of an oblate spheroid (oi = 02 ^ os) 

I-P 3 


Pi = P2 = 

P 3 = 


of 


(of - oi)3/2 [ 


(of — of)^/^ -|- Os arctan 


Os 


(of - oi)V2 


TT 

Os — 

2 


( 22 ) 

(23) 

(24) 

(25) 


Finally, in the case of a sphere (oi = 02 = os) one finish with Pi = P 2 = Ps = 4. 
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